An in-solution snapshot of SARS-COV-2 main protease maturation process and inhibition.

The main protease from SARS-CoV-2 (Mpro) is responsible for cleavage of the viral polyprotein. Mpro self-processing is called maturation, and it is crucial for enzyme dimerization and activity. Here we use C145S Mpro to study the structure and dynamics of N-terminal cleavage in solution. Native mass spectroscopy analysis shows that mixed oligomeric states are composed of cleaved and uncleaved particles, indicating that N-terminal processing is not critical for dimerization. A 3.5 Å cryo-EM structure provides details of Mpro N-terminal cleavage outside the constrains of crystal environment. We show that different classes of inhibitors shift the balance between oligomeric states. While non-covalent inhibitor MAT-POS-e194df51-1 prevents dimerization, the covalent inhibitor nirmatrelvir induces the conversion of monomers into dimers, even with intact N-termini. Our data indicates that the Mpro dimerization is triggered by induced fit due to covalent linkage during substrate processing rather than the N-terminal processing.

pH profiles of 3-chymotrypsin-like protease (3CLpro) from SARS-CoV-2 elucidate its catalytic mechanism and a histidine residue critical for activity.

3-Chymotrypsin-like protease (3CLpro) is a promising drug target for coronavirus disease 2019 and related coronavirus diseases because of the essential role of this protease in processing viral polyproteins after infection. Understanding the detailed catalytic mechanism of 3CLpro is essential for designing effective inhibitors of infection by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Molecular dynamics studies have suggested pH-dependent conformational changes of 3CLpro, but experimental pH profiles of SARS-CoV-2 3CLpro and analyses of the conserved active-site histidine residues have not been reported. In this work, pH-dependence studies of the kinetic parameters of SARS-CoV-2 3CLpro revealed a bell-shaped pH profile with 2 pKa values (6.9 ± 0.1 and 9.4 ± 0.1) attributable to ionization of the catalytic dyad His41 and Cys145, respectively. Our investigation of the roles of conserved active-site histidines showed that different amino acid substitutions of His163 produced inactive enzymes, indicating a key role of His163 in maintaining catalytically active SARS-CoV-2 3CLpro. By contrast, the H164A and H172A mutants retained 75% and 26% of the activity of WT, respectively. The alternative amino acid substitutions H172K and H172R did not recover the enzymatic activity, whereas H172Y restored activity to a level similar to that of the WT enzyme. The pH profiles of H164A, H172A, and H172Y were similar to those of the WT enzyme, with comparable pKa values for the catalytic dyad. Taken together, the experimental data support a general base mechanism of SARS-CoV-2 3CLpro and indicate that the neutral states of the catalytic dyad and active-site histidine residues are required for maximum enzyme activity.

Chemo- and Site-Selective Lysine Modification of Peptides and Proteins under Native Conditions Using the Water-Soluble Zolinium.

Site-selective lysine modification of peptides and proteins in aqueous solutions or in living cells is still a big challenge today. Here, we report a novel strategy to selectively quinolylate lysine residues of peptides and proteins under native conditions without any catalysts using our newly developed water-soluble zoliniums. The zoliniums could site-selectively quinolylate K350 of bovine serum albumin and inactivate SARS-CoV-2 3CLpro via covalently modifying two highly conserved lysine residues (K5 and K61). In living HepG2 cells, it was demonstrated that the simple zoliniums (5b and 5B) could quinolylate protein lysine residues mainly in the nucleus, cytosol, and cytoplasm, while the zolinium-fluorophore hybrid (8) showed specific lysosome-imaging ability. The specific chemoselectivity of the zoliniums for lysine was validated by a mixture of eight different amino acids, different peptides bearing potential reactive residues, and quantum chemistry calculations. This study offers a new way to design and develop lysine-targeted covalent ligands for specific application.

Discovery of 2-(furan-2-ylmethylene)hydrazine-1-carbothioamide derivatives as novel inhibitors of SARS-CoV-2 main protease.

The COVID-19 posed a serious threat to human life and health, and SARS-CoV-2 Mpro has been considered as an attractive drug target for the treatment of COVID-19. Herein, we report 2-(furan-2-ylmethylene)hydrazine-1-carbothioamide derivatives as novel inhibitors of SARS-CoV-2 Mpro developed by in-house library screening and biological evaluation. Similarity search led to the identification of compound F8-S43 with the enzymatic IC50 value of 10.76 µM. Further structure-based drug design and synthetic optimization uncovered compounds F8-B6 and F8-B22 as novel non-peptidomimetic inhibitors of Mpro with IC50 values of 1.57 µM and 1.55 µM, respectively. Moreover, enzymatic kinetic assay and mass spectrometry demonstrated that F8-B6 was a reversible covalent inhibitor of Mpro. Besides, F8-B6 showed low cytotoxicity with CC50 values of more than 100 µM in Vero and MDCK cells. Overall, these novel SARS-CoV-2 Mpro non-peptidomimetic inhibitors provide a useful starting point for further structural optimization.

Allosteric Hotspots in the Main Protease of SARS-CoV-2.

Inhibiting the main protease of SARS-CoV-2 is of great interest in tackling the COVID-19 pandemic caused by the virus. Most efforts have been centred on inhibiting the binding site of the enzyme. However, considering allosteric sites, distant from the active or orthosteric site, broadens the search space for drug candidates and confers the advantages of allosteric drug targeting. Here, we report the allosteric communication pathways in the main protease dimer by using two novel fully atomistic graph-theoretical methods: Bond-to-bond propensity, which has been previously successful in identifying allosteric sites in extensive benchmark data sets without a priori knowledge, and Markov transient analysis, which has previously aided in finding novel drug targets in catalytic protein families. Using statistical bootstrapping, we score the highest ranking sites against random sites at similar distances, and we identify four statistically significant putative allosteric sites as good candidates for alternative drug targeting.

Reconstruction of the unbinding pathways of noncovalent SARS-CoV and SARS-CoV-2 3CLpro inhibitors using unbiased molecular dynamics simulations.

The main protease (3CLpro) is one of the essential components of the SARS-CoVs viral life cycle, which makes it an interesting target for overpowering these viruses. Although many covalent and noncovalent inhibitors have been designed to inhibit this molecular target, none have gained FDA approval as a drug. Because of the high rate of COVID-19 pandemic development, in addition to laboratory research, we require in silico methods to accelerate rational drug design. The unbinding pathways of two SARS-CoV and SARS-CoV-2 3CLpro noncovalent inhibitors with the PDB IDs: 3V3M, 4MDS, 6W63, 5RF7 were explored from a comparative perspective using unbiased molecular dynamics (UMD) simulations. We uncovered common weak points for selected inhibitors that could not interact significantly with a binding pocket at specific residues by all their fragments. So water molecules entered the free binding S regions and weakened protein-inhibitor fundamental interactions gradually. N142, G143, and H163 are the essential residues, which cause key protein-ligand interactions in the binding pocket. We believe that these results will help design new potent inhibitors against SARS-CoV-2.

Surface cysteines could protect the SARS-CoV-2 main protease from oxidative damage.

The SARS-CoV-2 main protease (Mpro) is responsible for cleaving twelve nonstructural proteins from the viral polyprotein. Mpro, a cysteine protease, is characterized by a large number of noncatalytic cysteine (Cys) residues, none involved in disulfide bonds. In the absence of a tertiary-structure stabilizing role for these residues, a possible alternative is that they are involved in redox processes. We report experimental work in support of a proposal that surface cysteines on Mpro can protect the active-site Cys145 from oxidation by reactive oxygen species (ROS). In investigations of enzyme kinetics, we found that mutating three surface cysteines to serines did not greatly affect activity, which in turn indicates that these cysteines could protect Cys145 from oxidative damage.

Optimising Elastic Network Models for Protein Dynamics and Allostery: Spatial and Modal Cut-offs and Backbone Stiffness.

The family of coarse-grained models for protein dynamics known as Elastic Network Models (ENMs) require careful choice of parameters to represent well experimental measurements or fully-atomistic simulations. The most basic ENM that represents each protein residue by a node at the position of its C-alpha atom, all connected by springs of equal stiffness, up to a cut-off in distance. Even at this level a choice is required of the optimum cut-off distance and the upper limit of elastic normal modes taken in any sum for physical properties, such as dynamic correlation or allosteric effects on binding. Additionally, backbone-enhanced ENM (BENM) may improve the model by allocating a higher stiffness to springs that connect along the protein backbone. This work reports on the effect of varying these three parameters (distance and mode cutoffs, backbone stiffness) on the dynamical structure of three proteins, Catabolite Activator Protein (CAP), Glutathione S-transferase (GST), and the SARS-CoV-2 Main Protease (M pro ). Our main results are: (1) balancing B-factor and dispersion-relation predictions, a near-universal optimal value of 8.5 Å is advisable for ENMs; (2) inhomogeneity in elasticity brings the first mode containing spatial structure not well-resolved by the ENM typically within the first 20; (3) the BENM only affects modes in the upper third of the distribution, and, additionally to the ENM, is only able to model the dispersion curve better in this vicinity; (4) BENM does not typically affect fluctuation-allostery, which also requires careful treatment of the effector binding to the host protein to capture.

Structural Basis of Main Proteases of Coronavirus Bound to Drug Candidate PF-07304814.

New variants of the severe acute respiratory syndrome Coronavirus 2 (SARS-CoV-2) emerged and spread rapidly all over the world, which strongly supports the need for pharmacological options to complement vaccine strategies. Main protease (Mpro or 3CLpro) is a critical enzyme in the life cycle of SARS-CoV-2 and appears to be highly conserved among different genera of coronaviruses, making it an ideal target for the development of drugs with broad-spectrum property. PF-07304814 developed by Pfizer is an intravenously administered inhibitor targeting SARS-CoV-2 Mpro. Here we showed that PF-07304814 displays broad-spectrum inhibitory activity against Mpros from multiple coronaviruses. Crystal structures of Mpros of SARS-CoV-2, SARS-CoV, MERS-CoV, and HCoV-NL63 bound to the inhibitor PF-07304814 revealed a conserved ligand-binding site, providing new insights into the mechanism of inhibition of viral replication. A detailed analysis of these crystal structures complemented by comprehensive comparison defined the key structural determinants essential for inhibition and illustrated the binding mode of action of Mpros from different coronaviruses. In view of the importance of Mpro for the medications of SARS-CoV-2 infection, insights derived from the present study should accelerate the design of pan-coronaviral main protease inhibitors that are safer and more effective.

Effect of pH on stability of dimer structure of the main protease of coronavirus-2.

The viral main protease (Mpro) from a novel severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) is a key enzyme essential for viral replication and has become an attractive target for antiviral drug development. The Mpro forms a functional dimer and exhibits a pH-dependent enzyme activity and dimerization. Here, we report a molecular dynamics (MD) investigation to gain insights into the structural stability of the enzyme dimer at neutral and acidic pH. Our data shows larger changes in structure of the protein with the acidic pH than that with the neutral pH. Structural analysis of MD trajectories reveals a substantial increase in intersubunit separation, the loss of domain contacts, binding free energy and interaction energy of the dimer which implies the protein instability and tendency of dimer dissociation at acidic pH. The loss in the interaction energy is mainly driven by electrostatic interactions. We have identified the intersubunit hydrogen-bonding residues involved in the decreased dimer stability. These findings may be helpful for rational drug design and target evaluation against COVID-19.

Structural basis for replicase polyprotein cleavage and substrate specificity of main protease from SARS-CoV-2.

The main protease (Mpro) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a key enzyme, which extensively digests CoV replicase polyproteins essential for viral replication and transcription, making it an attractive target for antiviral drug development. However, the molecular mechanism of how Mpro of SARS-CoV-2 digests replicase polyproteins, releasing the nonstructural proteins (nsps), and its substrate specificity remain largely unknown. Here, we determine the high-resolution structures of SARS-CoV-2 Mpro in its resting state, precleavage state, and postcleavage state, constituting a full cycle of substrate cleavage. The structures show the delicate conformational changes that occur during polyprotein processing. Further, we solve the structures of the SARS-CoV-2 Mpro mutant (H41A) in complex with six native cleavage substrates from replicase polyproteins, and demonstrate that SARS-CoV-2 Mpro can recognize sequences as long as 10 residues but only have special selectivity for four subsites. These structural data provide a basis to develop potent new inhibitors against SARS-CoV-2.

A new inactive conformation of SARS-CoV-2 main protease.

The SARS-CoV-2 main protease (Mpro) has a pivotal role in mediating viral genome replication and transcription of the coronavirus, making it a promising target for drugs against the COVID-19 pandemic. Here, a crystal structure is presented in which Mpro adopts an inactive state that has never been observed before, called new-inactive. It is shown that the oxyanion loop, which is involved in substrate recognition and enzymatic activity, adopts a new catalytically incompetent conformation and that many of the key interactions of the active conformation of the enzyme around the active site are lost. Solvation/desolvation energetic contributions play an important role in the transition from the inactive to the active state, with Phe140 moving from an exposed to a buried environment and Asn142 moving from a buried environment to an exposed environment. In new-inactive Mpro a new cavity is present near the S2' subsite, and the N-terminal and C-terminal tails, as well as the dimeric interface, are perturbed, with partial destabilization of the dimeric assembly. This novel conformation is relevant both for comprehension of the mechanism of action of Mpro within the catalytic cycle and for the successful structure-based drug design of antiviral drugs.

An extended conformation of SARS-CoV-2 main protease reveals allosteric targets.

SignificanceThe coronavirus main protease (Mpro) is required for viral replication. Here, we obtained the extended conformation of the native monomer of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Mpro by trapping it with nanobodies and found that the catalytic domain and the helix domain dissociate, revealing allosteric targets. Another monomeric state is termed compact conformation and is similar to one protomer of the dimeric form. We designed a Nanoluc Binary Techonology (NanoBiT)-based high-throughput allosteric inhibitor assay based on structural conformational change. Our results provide insight into the maturation, dimerization, and catalysis of the coronavirus Mpro and pave a way to develop an anticoronaviral drug through targeting the maturation process to inhibit the autocleavage of Mpro.

Modulation of the monomer-dimer equilibrium and catalytic activity of SARS-CoV-2 main protease by a transition-state analog inhibitor.

The role of dimer formation for the onset of catalytic activity of SARS-CoV-2 main protease (MProWT) was assessed using a predominantly monomeric mutant (MProM). Rates of MProWT and MProM catalyzed hydrolyses display substrate saturation kinetics and second-order dependency on the protein concentration. The addition of the prodrug GC376, an inhibitor of MProWT, to MProM leads to an increase in the dimer population and catalytic activity with increasing inhibitor concentration. The activity reaches a maximum corresponding to a dimer population in which one active site is occupied by the inhibitor and the other is available for catalytic activity. This phase is followed by a decrease in catalytic activity due to the inhibitor competing with the substrate. Detailed kinetics and equilibrium analyses are presented and a modified Michaelis-Menten equation accounts for the results. These observations provide conclusive evidence that dimer formation is coupled to catalytic activity represented by two equivalent active sites.

Improved SARS-CoV-2 main protease high-throughput screening assay using a 5-carboxyfluorescein substrate.

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) as a global threat to human health has highlighted the need for the development of novel therapies targeting current and emerging coronaviruses with pandemic potential. The coronavirus main protease (Mpro, also called 3CLpro) is a validated drug target against coronaviruses and has been heavily studied since the emergence of SARS-CoV-2 in late 2019. Here, we report the biophysical and enzymatic characterization of native Mpro, then characterize the steady-state kinetics of several commonly used FRET substrates, fluorogenic substrates, and six of the 11 reported SARS-CoV-2 polyprotein cleavage sequences. We then assessed the suitability of these substrates for high-throughput screening. Guided by our assessment of these substrates, we developed an improved 5-carboxyfluorescein-based FRET substrate, which is better suited for high-throughput screening and is less susceptible to interference and false positives than existing substrates. This study provides a useful framework for the design of coronavirus Mpro enzyme assays to facilitate the discovery and development of therapies targeting Mpro.

Ultralarge Virtual Screening Identifies SARS-CoV-2 Main Protease Inhibitors with Broad-Spectrum Activity against Coronaviruses.

Drugs targeting SARS-CoV-2 could have saved millions of lives during the COVID-19 pandemic, and it is now crucial to develop inhibitors of coronavirus replication in preparation for future outbreaks. We explored two virtual screening strategies to find inhibitors of the SARS-CoV-2 main protease in ultralarge chemical libraries. First, structure-based docking was used to screen a diverse library of 235 million virtual compounds against the active site. One hundred top-ranked compounds were tested in binding and enzymatic assays. Second, a fragment discovered by crystallographic screening was optimized guided by docking of millions of elaborated molecules and experimental testing of 93 compounds. Three inhibitors were identified in the first library screen, and five of the selected fragment elaborations showed inhibitory effects. Crystal structures of target-inhibitor complexes confirmed docking predictions and guided hit-to-lead optimization, resulting in a noncovalent main protease inhibitor with nanomolar affinity, a promising in vitro pharmacokinetic profile, and broad-spectrum antiviral effect in infected cells.

Mechanism of Caspase-1 Inhibition by Four Anti-inflammatory Drugs Used in COVID-19 Treatment.

The inflammatory protease caspase-1 is associated with the release of cytokines. An excessive number of cytokines (a "cytokine storm") is a dangerous consequence of COVID-19 infection and has been indicated as being among the causes of death by COVID-19. The anti-inflammatory drug colchicine (which is reported in the literature to be a caspase-1 inhibitor) and the corticosteroid drugs, dexamethasone and methylprednisolone, are among the most effective active compounds for COVID-19 treatment. The SERM raloxifene has also been used as a repurposed drug in COVID-19 therapy. In this study, inhibition of caspase-1 by these four compounds was analyzed using computational methods. Our aim was to see if the inhibition of caspase-1, an important biomolecule in the inflammatory response that triggers cytokine release, could shed light on how these drugs help to alleviate excessive cytokine production. We also measured the antioxidant activities of dexamethasone and colchicine when scavenging the superoxide radical using cyclic voltammetry methods. The experimental findings are associated with caspase-1 active site affinity towards these compounds. In evaluating our computational and experimental results, we here formulate a mechanism for caspase-1 inhibition by these drugs, which involves the active site amino acid Cys285 residue and is mediated by a transfer of protons, involving His237 and Ser339. It is proposed that the molecular moiety targeted by all of these drugs is a carbonyl group which establishes a S(Cys285)-C(carbonyl) covalent bond.

Structural and energetic features of the dimerization of the main proteinase of SARS-CoV-2 using molecular dynamic simulations.

The COVID-19 pandemic caused by SARS-CoV-2 has been declared a global health crisis. The development of anti-SARS-CoV-2 drugs heavily depends on the systematic study of the critical biological processes of key proteins of coronavirus among which the main proteinase (Mpro) dimerization is a key step for virus maturation. Because inhibiting the Mpro dimerization can efficiently suppress virus maturation, the key residues that mediate dimerization can be treated as targets of drug and antibody developments. In this work, the structure and energy features of the Mpro dimer of SARS-CoV-2 and SARS-CoV were studied using molecular dynamics (MD) simulations. The free energy calculations using the Generalized Born (GB) model showed that the dimerization free energy of the SARS-CoV-2 Mpro dimer (-107.5 ± 10.89 kcal mol-1) is larger than that of the SARS-CoV Mpro dimer (-92.83 ± 9.81 kcal mol-1), indicating a more stable and possibly a quicker formation of the Mpro dimer of SARS-CoV-2. In addition, the energy decomposition of each residue revealed 11 key attractive residues. Furthermore, Thr285Ala weakens the steric hindrance between the two protomers of SARS-CoV-2 that can form more intimate interactions. It is interesting to find 11 repulsive residues which effectively inhibit the dimerization process. At the interface of the Mpro dimer, we detected three regions that are rich in interfacial water which stabilize the SARS-CoV-2 Mpro dimer by forming hydrogen bonds with two protomers. The key residues and rich water regions provide important targets for the future design of anti-SARS-CoV-2 drugs through inhibiting Mpro dimerization.

Structural insights into the substrate-binding site of main protease for the structure-based COVID-19 drug discovery.

An attractive drug target to combat COVID-19 is the main protease (Mpro ) because of its key role in the viral life cycle by processing the polyproteins translated from the viral RNA. Studying the crystal structures of the protease is important to enhance our understanding of its mechanism of action at the atomic-level resolution, and consequently may provide crucial structural insights for structure-based drug discovery. In the current study, we report a comparative structural analysis of the Mpro substrate binding site for both apo and holo forms to identify key interacting residues and conserved water molecules during the ligand-binding process. It is shown that in addition to the catalytic dyad residues (His41 and Cys145), the oxyanion hole residues (Asn142-Ser144) and residues His164-Glu166 form essential parts of the substrate-binding pocket of the protease in the binding process. Furthermore, we address the issue of the substrate-binding pocket flexibility and show that two adjacent loops in the Mpro structures (residues Thr45-Met49 and Arg188-Ala191) with high flexibility can regulate the binding cavity' accessibility for different ligand sizes. Moreover, we discuss in detail the various structural and functional roles of several important conserved and mobile water molecules within and around the binding site in the proper enzymatic function of Mpro . We also present a new docking protocol in the framework of the ensemble docking strategy. The performance of the docking protocol has been evaluated in predicting ligand binding pose and affinity ranking for two popular docking programs; AutoDock4 and AutoDock Vina. Our docking results suggest that the top-ranked poses of the most populated clusters obtained by AutoDock Vina are the most important representative docking runs that show a very good performance in estimating experimental binding poses and affinity ranking.

Protocol for structure determination of SARS-CoV-2 main protease at near-physiological-temperature by serial femtosecond crystallography.

The SARS-CoV-2 main protease of (Mpro) is an important target for SARS-CoV-2 related drug repurposing and development studies. Here, we describe the steps for structural characterization of SARS-CoV-2 Mpro, starting from plasmid preparation and protein purification. We detail the steps for crystallization using the sitting drop, microbatch (under oil) approach. Finally, we cover data collection and structure determination using serial femtosecond crystallography. For complete details on the use and execution of this protocol, please refer to Durdagi et al. (2021).